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Abstract 

We consider the discrete solitons bifurcating from the anti-continuum limit of the discrete 
nonlinear Schrodinger (NLS) lattice. The discrete soliton in the anti-continuum limit repre- 
sents an arbitrary finite superposition of in-phase or anti-phase excited nodes, separated by 
an arbitrary sequence of empty nodes. By using stability analysis, we prove that the discrete 
solitons are all unstable near the anti-continuum limit, except for the solitons, which consist 
of alternating anti-phase excited nodes. We classify analytically and confirm numerically the 
number of unstable eigenvalues associated with each family of the discrete solitons. 



1 Introduction 

Nonlinear instabilities and emergence of coherent structures in difFerential-difFerence equations have 
become topics of physical importance and mathematical interest in the past decade. Numerous 
applications of these problems have emerged ranging from nonlinear optics, in the dynamics of 
guided waves in inhomogeneous optical structures ^1 H] and photonic crystal lattices O 0], to 
atomic physics, in the dynamics of Bose-Einstein condensate droplets in periodic (optical lattice) 
potentials [HIEIIZIIHI and from condensed matter, in Josephson-junction ladders OEIl, to biophysics, 
m various models of the DNA double strand [TTlini- This large range of models and applications 
has been summarized by now in a variety of reviews such as |13l 1141 1151 1161 117j . 

One of the prototypical differential-difference models that is both physically relevant and mathe- 
matically tractable is the so-called discrete nonlinear Schrodinger (NLS) equation, 

iUn+ (3/^dUn+l\Un?Un^Q, (1.1) 

where u„ = Un{t) is a complex amplitude in time i, n G Z'' is the d-dimensional lattice, is the 
d-dimensional discrete Laplacian, (3 is the dispersion coefficient, and 7 is the nonlinearity coefficient. 
Before we delve into mathematical analysis of the discrete NLS equation, it would be relevant to 
discuss briefly the recent physical applications of this model. 

The most direct implementation of the discrete NLS equation can be identified in one-dimensional 
arrays of coupled optical waveguides |5j . These may be multi-core structures created in a slab 
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of a semiconductor material (such as AlGaAs), or virtual ones, induced by a set of laser beams 
illuminating a photorefractive crystal. In this experimental implementation, there are about forty 
lattice sites (guiding cores), and the localized modes (discrete solitons) may propagate over twenty 
diffraction lengths. 

Light-induced photonic lattices E] have recently emerged as another application of the discrete 
NLS equation. The refractive index of the nonlinear photonic lattices changes periodically due to a 
grid of strong beams, while a weaker probe beam is used to monitor the localized modes (discrete 
solitons). A number of promising experimental studies of discrete solitons in light-induced photonic 
lattices was reported recently in physics literature. 

An array of Bose-Einstein condensate droplets trapped in a strong optical lattice with thousands 

of atoms in each droplet, is another direct physical realization of the discrete NLS equation 

In this context, the model can be derived systematically by using the Wannier function expansions 

13 ig. 

Besides applications to optical waveguides, photonic crystal lattices, and Bose-Einstein condensates 
trapped in optical lattices, the discrete NLS equation also arises as the envelope wave reduction of 
the general nonlinear Klein-Gordon lattices 118'. 

This rich variety of physical contexts makes it timely and relevant to analyze the mathematical 
aspects of the discrete NLS equation 1)1.1(1 , including the existence and stability of localized modes 
(discrete solitons). A very helpful tool for such an analysis is the so-called anti-continuum limit 
P ^ ^3 1 where the nonlinear oscillators of the model are uncoupled. Existence of localized 
modes in this limit can be characterized in full details jl5| . Persistence, multiplicity, and stability 
of localized modes can be studied with continuation methods both analytically and numerically 

m 

We study localized modes of the discrete NLS equation 1(1. l|l in the forthcoming scries of two papers. 
This first paper describes stability analysis of discrete solitons in the onc-dimcnsional NLS lattice 
{d = 1). The second paper will present Lyapunov-Schmidt reductions for persistence, multiplicity 
and stability of discrete vortices in the two-dimensional NLS lattice {d — 2). 

This paper is structured as follows. We review the known results on existence of one-dimensional 
discrete solitons in Section 2. General stability and instability results for discrete solitons in the 
anti-continuum limit are proved in Section 3. These stability results are illustrated for two particular 
families of the discrete solitons in Sections 4 and 5. Besides explicit perturbation series expansions 
results, we compare asymptotic approximations and numerical computations of stable and unstable 
eigenvalues in the linearized stability problem. Section 6 concludes the first paper of this series. 

2 Existence of discrete solitons 

We consider the normalized form of the discrete NLS equation 1(1.1(1 in one dimension {d — 1): 

iiln + e (u„+l - 2Un + Un-l) + |M,ipu„ = 0, (2.1) 

where w„(i) : M-|_ ^ C, n G Z, and e > is the inverse squared step size of the discrete one- 
dimensional NLS lattice. The discrete solitons are given by the time-periodic solutions of the 
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discrete NLS equation H2.1|l : 

-0„e^^^"''^*+*'", MeM, 0„eC, neZ, (2.2) 
where G is parameter and (/i, (/)„) solve the nonhnear difference equations on ri G Z: 

(^ - |0„n0n = e (0n+l + (l)n-l) ■ (2.3) 

Existence of the discrete solitons was studied recently in j20[l21[|^ . inspired by the pioneer papers 
|23l . Recent summary of the existence results is given in JHI ■ Since discrete solitons in the 
focusing NLS lattice 1)2. l|l may exist only for fi > 19 and the parameter fi is scaled out by the 
scaling transformation, 

0« = VT'^n, £ = (2.4) 

the parameter fj, > will henceforth be set to fi — I. Another arbitrary parameter 9o, which exists 
due to the gauge invariance of the discrete NLS equation 1)2.111 , is incorporated in the anzats 1)2.2(1 , 
such that at least one value of 0„ can be chosen real-valued without lack of generality. Using this 
convention, we represent below the known existence results. 

Proposition 2.1 There exist eo > 0, k > and 4>ao > 0, such that the difference equations ^2.S\) 
with jJ, — I and < e < cq have continuous families of the discrete solitons with the properties: 

" ^' (2.5) 

I e Z\S, ^ ' 

(2.6) 

(2.7) 

where S is a finite set of nodes of the lattice n G Z and 0„ — {0, tt}, n £ S . 

Proof. See Theorem 2.1 and Appendices A and B in |12| for the proof of the limiting solution 
()2.5)l from the inverse function theorem. See Theorem 3 in |23 for the proof of the exponential 
decay ()2.6)) from the bound estimates. See Section 3.2 in |^ for the proof of the reality condition 
()2.7)l from the conservation of the density current. Various theoretical and numerical bounds on eo 
are obtained in [13 [13 EDI El] . ■ 

Due to the property 1)2. 7)1 . the difference equations 1)2.3)1 with /i = 1 can be rewritten as follows: 

{l-<j,l)<j>n^e{<l>„+i+<j,„^i), neZ. (2.8) 
For our analysis, we shall derive two technical results on properties of solutions (/)„, n € Z. 

Lemma 2.2 There exists eo > 0, such that the solution n € Z is represented by the convergent 
power series for < e < cq: 

OG 

k=l 

where (j)n^ is given by 1^2. 5\l . 
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Proof. The statement follows from the Implicit Function Theorem (see Theorem 2.7.2 in 



since the Jacobian matrix for the system H2.8|l is non-singular at (/>„ 
right-hand-side of the system H2.8|l is analytic in e. 



, n e Z, while the 



Lemma 2.3 There exists < ei < ep, such that the number of changes in the sign oj <j)n on n ^"L 
for < e < ei equals to the number of ti- differences of the adjacent On, n S in the limiting 
solution 1^2. 5\} . 

Proof. Consider two adjacent excited nodes ni,n2 € S, separated by N empty nodes, such that 
712 — ni — 1 + N and iV > 1. We need to prove that the number of 7r-differences in the argument 
(j)n, ni < n < n2 for small e > is exactly one if — 9^ — tt and zero if ~ ^ni — 0- To do 
so, we consider the difference equations (12. 8|) on ni < n < 712 as the N-hy-N matrix system: 



N 



eb 



(2.10) 



where 



N 



A 



N 



I 1 



V 







-1)-'", and 











1 



3 AT 



-1 / 





V <^n2 } 



(2.11) 



Let -D/.j, l</<J<^be the determinant of the block of the matrix yljv between /-th and J-th 
rows and columns. By Cramer's rule, we have 



N 



(2.12) 



Since limj^o Di j = 1 for all 1 < / < J < A^, we have 



lim e 

lim e^-'f 
lim e-''-i-^( 



'ni+j 



ni+j 



1 < J < 



1 

N + 1 



On2, J = 

N + 1 



2 ' 
<j<N. 



The statement of Lemma follows from the signs of 0„, ni < n < n2 for small e > 0. 



By Proposition and Lemma all families of the discrete solitons as e ^ can be classified by 
a sequence of {0}, {-I-}, and {— } of the limiting solution H2.F)|1 on the finite set 5* In particular, 
we consider two ordered sets S: 

5i = {1,2,3, ...,iV} (2.13) 

and 

^2 = {1,3,5,. ..,2iV-l}, (2.14) 

where dim(S'i) = dim(52) = iV < 00. The set Si includes the Page mode {N = 2: 9i = 82 = 0) 
and the twisted mode {N — 2: 9i ~ 0, ^2 ~ tt). The set S2 includes the Page and twisted modes 
{N = 2), separated by an empty node. 
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3 Stability of discrete solitons 



The spectral stability of discrete solitons is studied with the standard linearization: 

w„W = 6^^'"''^*+'''' (0n + a«e^*+5„e^*) , A G C, (a„,5„)eC2, n e Z, (3.1) 

where (A,a„,5„) solve the linear eigenvalue problem on n G Z: 

(l ~ 20^) a„ - (?!)^6„ - e (a„+i + a„_i) = iAa„, 
-<l>lan + (1 - 2(1)1) bn - e {bn+i + = -iXbn- (3.2) 

The discrete soliton (|2.2|l is called spectrally unstable if there exists A and {an,bn), n G Z in the 
problem H3.2() . such that Re(A) > and J2nez (knP + l^nP) < c)o. Otherwise, the soliton is called 
weakly spectrally stable. Orbital stability of the discrete one-pulse soliton was studied in the anti- 
continuum limit e ^ |2f)j and close to the contimmm limit e ^ oo 27 . Spectral instabilities of two- 
pulse and multi-pulse solitons were considered in |28[ I29L ISUl 1311 IS^ I33j by numerical and variational 
approximations. It was well understood from intuition supported by numerical simulations )13l I32| 
that the discrete solitons with the alternating sequence of 0„ = {0, tt} in the limiting solution H2.5|l 
are spectrally stable as e — > but have eigenvalues with so-called negative Krein signature, which 
become complex by means of the Hamiltonian-Hopf bifurcations [211 EHl- AH other families of 
discrete solitons have unstable real eigenvalues A in the anti-continuum limit for any e 7^ 

Here we prove these preliminary observations and find the precise number of stable and unstable 
eigenvalues in the linearized stability problem H3.2|l for small e > 0. Our results are similar to 
those in the Lyapunov-Schmidt reductions, which are applied to continuous multi-pulse solitons in 
nonlinear Schrodinger equations |34l 1351 136[ I37| . In particular, the main conclusion on stability 
of alternating up-down solitons and instability of any other up-up and down-down sequences of 
solitons was found for multi-pulse homoclinic orbits arising in the so-called orbit-flip bifurcation 
|34l p. 176]. The same conclusion agrees with qualitative predictions for the discrete NLS equations 
p.66]. 

Let fl — L^(Z, C) be the Hilbert space of square-summable bi-infinite complex- valued sequences 
{u„}„gz, equipped with the inner product and norm: 

(u,w)o ^u„w„, ||u||^ = ^ |u„p <- (-3 3) 

We use bold notations u for an infinite-dimensional vector in that consists of components Un for 
all rt e Z. The stability problem (|3.2|) is transformed with the substitution, 

an^Un + iWn, b^ ^ - iw^, n e Z, (3.4) 

to the form: 

(1 - 3(/)^) u„ - e (u„+i + = -Au;„, 

{1 - (1)1) Wn - e [Wn+l + Wn~l) = Am„. (3.5) 

The matrix- vector form of the problem (|3.5|l is 

£+u = —Aw, C-w = Au, (3.6) 
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where C± are infinite-dimensional symmetric tri-diagonal matrices, which consist of elements: 



('C+)„,„ = 1 - 30„, (£_)„„ = i-0„, (/:±)„^„+i = (/:±)„+i ,„ = 

Equivalently, the stability problem H3.6() is rewritten in the Hamiltonian form: 

JHil^ ^ Xip, (3.7) 

where ^p is the infinite-dimensional eigenvector, which consists of 2-blocks of (u„,w„)"^, is the 
infinite-dimensional skew-symmetric matrix, which consists of 2-by-2 blocks of 




and Ti. is the infinite-dimensional symmetric matrix, which consists of 2-by-2 blocks of 



('^+)n.m 





The representation H3.7|l follows from the Hamiltonian structure of the discrete NLS equation (|2.1|l . 
where J7 is the symplectic operator and Ti. is the linearized Hamiltonian. By Lemma l2.2l the matrix 
Ti. is expanded into the power series: 




H = -f ^ e*'-H('=\ (3.8) 
fe=i 

where Ti^^^ is diagonal with two blocks: 

/ 1 n \ 

(3.9) 

Let N ~ dim(S') < oo. The spectrum of Ti'^'^^ip = 7¥> has exactly N negative eigenvalues 7 — —2, 
N zero eigenvalues 7 = and infinitely many positive eigenvalues 7 = +1. The negative and zero 
eigenvalues 7 = —2 and 7 = map to N double zero eigenvalues A = in the eigenvalue problem 
JTi'^^^xj} = Xij}. The positive eigenvalues 7 = -f 1 map to the infinitely many eigenvalues A — ±i. 

Since finitely many zero eigenvalues of JTi^^^ are isolated from the rest of the spectrum, their shifts 
vanish as e — > 0, according to the regular perturbation theory '^S'. We can therefore locate small 
unstable eigenvalues Re(A) > of the stability problem (|3.7|) for small e > from their limits at 
6 = 0. On the other hand, infinitely many imaginary eigenvalues of j/Ti'"' become the continuous 
spectrum band as e 7^ |39| . However, since the difference operator J^Ti has exponentially decaying 
potentials (j>n, n G Z, due to the decay condition (|2.6|l . the continuous spectral bands of J^Ti are 
located on the imaginary axis of A near the points A = ±i, similarly to the case = 0, n € Z |39| . 
Therefore, the infinite-dimensional part of the spectrum does not produce any unstable eigenvalues 
Re(A) > in the stability problem 1)3. 7|l as e > 0. Results of the regular perturbation theory are 
formulated and proved below. 

Lemma 3.1 Assume that n ^ Z is the discrete soliton, described in Provosition \2.1\ Let 
N — dim(S') < 00. Let jj, 1 < j < N be small eigenvalues of Ti as e — > 0, such that 

lini7j=0, l<i<A^. (3.10) 
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There exists < e* < eo, such that the eigenvalue problem jcl.T)) with n S Z and < e < has 
N pairs of small eigenvalues Xj and —Xj, ^ l£ j l£ N , that satisfy the leading-order behavior: 

lim ^ = 2, 1< 7 < N. (3.11) 

Proof. Since the operator is Fredholm of zero index and empty kernel at e = 0, it can be 
inverted for small e > and the non-self-adjoint eigenvalue problem (|3.6|l can be transformed to 
the self-adjoint diagonalization problem: 

/:_w = -A2/:;1w, (3.12) 

such that 

- -P^^^, (3.13) 
(w,£_|_ w)o 

where the inner product is defined in l|3.3|l . Since all small eigenvalues of Ti. are small eigenvalues of 
we denote Wj be an eigenvector of which corresponds to the small eigenvalue , I < j < N 
in the limiting condition H3.10|l . By continuity of the eigenvectors and completeness of kcr(£'^-'), 
there exists a set of normalized coefficients {cn.j}nes for each 1 < j < N, such that 

lim wj = wj"' = c„je„, |c„jp = 1, (3.14) 



where e„ is the unit vector in Q. It follows from the direct computations that 

lim(w„/:;iw,) = (wf\4")-Vf)) = -i. (3.15) 



The leading-order behavior (|3.11() follows from H3.13|l and H3.15|l by the regular perturbation theory 

m ■ 

Corollary 3.2 Each small positive eigenvalue "fj corresponds to a pair of positive and negative 
eigenvalues Xj and ~Xj for small e > 0. Each small negative eigenvalue "/j corresponds to a pair 
of purely imaginary eigenvalues Xj and —Xj for small e > 0. The latter eigenvalues have negative 
Krein signature: 

Hip) = (u, C+u) + (w, £_ w) = 2 (w, /:_w) < 0. (3.16) 

For any e 7^ 0, there exists a simple zero eigenvalue of Ti. due to the gauge symmetry of the discrete 
solitons H2.2|l . as the parameter Oq is arbitrary, such that — 0. When all other {N — 1) 

eigenvalues 7j are non-zero for any e 7^ 0, the splitting of the semi-simple zero eigenvalue of Ti^"-* 
is called generic. The generic splitting gives a sufficient condition for unique (up to the gauge 
invariance) continuation of discrete solitons for e [H^j which is also garanteed by Proposition 

OEl 

Let no and po be the numbers of negative and positive eigenvalues 7^, defined in Lemma l3. II The 
splitting is generic if po = — 1 — no. The numbers no and po are computed exactly from the 
limiting solution (|2.5|l as follows. 

Lemma 3.3 There exists < ei < eo; such that the index no for < e < ei equals to the number 
of TT- differences of the adjacent On, n E S in the limiting solution ^g..5|) . while po = N — 1 — uq. 



7 



Proof. Since = for any < e < eo, the number no of negative eigenvalues of £_ coincides 
with the number of times when <p changes the sign, by the Discrete Sturm-Liouville Theorem |39j . 
In the case e = 0, this number equals the number of 7r-dilferences of the adjacent 6'„, n G S in the 
limiting solution (|2.5I) . By Lemma |2.3I the number remains continuous as e 7^ 0. The difference 
equation >C_w = has only two fundamental solutions, such that w = ciWi + C2W2, where Ci,C2 
are arbitrary parameters, wi = </> is exponentially decaying as \n\ — > 00, and W2 is exponentially 
growing as |n| 00, due to the discrete Wronskian identity As a result, the kernel of £_ is 
one-dimensional for e 7^ 0, such that po — N — 1 — no- ■ 

It was recently studied 001 ^| that there exists a closure relation between the negative index of the 
linearized Hamiltonian Ti. and the number of unstable eigenvalues of the linearized operator J^Ti. 
The closure relation can be extended from the coupled NLS equations to the discrete NLS equations 
by using the same methods |40[ HT] . We hence formulate the closure relation for the discrete NLS 
equations p. 1(1 . 

Proposition 3.4 Let niTL) he the finite number of negative eigenvalues of Ti.. Let iVroai be the 
number of positive real eigenvalues A in the problem jS.T)) , N^^^ be the number of pairs of purely 
imaginary eigenvalues A with negative Krein signature {ip,TL'ip) < 0, and A^comp be the number of 
complex eigenvalues A in the first open quadrant of A. Let p{P') = 1 if P' > and p{P') = if 
P' < 0, where 

P' = ml-e-^jml (3.17) 

Assume that X ~ is a double eigenvalue of the problem 7| j. Assume that no purely imaginary 
eigenvalues A exist inside the continuous spectrum or have zero Krein signature. The indices above 
satisfy the closure relation: 

n{n) - p{P') = iVreal + 2Ai- + 2iV,o„,p. (3.18) 

Proof. The left-hand-side of H3.18|l is the negative index of Ti in the constrained subspace of 
J7, which is reduced by one, if the power \\(f>\\Q is increasing function of fi. Due to the scaling 
transformation (|2.4() . the derivative of ||0||q in /i is given by H3.17|l . where the hats for and e are 
omitted. The right-hand-side of (|3.18|) is the negative index of H on the subspace of il, associated to 
the eigenvalue problem (|3.t)|) . The two indices are equal under the assumptions of the proposition, 
according to gOl ^ • ■ 

Corollary 3.5 There exists < 62 < ei, such that the indices of Proposition for < e < 62 
equal to the indices of Lemmas \3.1\ and \S.3[ as follows: 

n{n)=N + no, p{P') = l, iV,eal = A - 1 - no, A^i^nag = '^O, ^comp=0, (3.19) 

and the closure relation J^S.l^) is met. 

When the assumptions of Proposition 13.41 are not satisfied, instability bifurcations may occur in 
the eigenvalue problem H3.7|l . which results in the redistribution of the numbers n(7i), p{P'), A'lcai, 
-^imagi Acomp- The Hamiltouian-Hopf bifurcation, which is typical for the discrete multi- 
humped solitons [23 EH El, occurs when the purely imaginary eigenvalues A of negative Krein 
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signature (■0,7^-0) < collide with the purely imaginary eigenvalues A of positive Krein signature 
{tpjT-lxp) > or with the continuous spectral band and bifurcate as complex unstable eigenvalues A 
with Re(A) > 0. It follows from Corollaries 13 . 21 and 13 . 51 that there can be at most uq Hamiltonian- 
Hopf instability bifurcations, which result in at most N + tiq — 1 unstable eigenvalues, unless the 
indices n(7i) and p{P') change as a result of the zero eigenvalue bifurcations. 

Combining Lemmas 13.11 and 13.31 and Corollaries 13.21 and 13.51 we summarize the main stability- 
instability result for the discrete solitons of the discrete NLS equation (|2.1|l . 

Theorem 3.6 Let uq be the number of n- differences of the adjacent 6n, n S in the limiting 
solution 1^2. 5\l . The discrete soliton is spectrally stable for small e > if and only if hq ^ N — 
1. When uq < N — 1, the discrete soliton is spectrally unstable with exactly N — 1 — uq real 
unstable eigenvalues A in the problem jS. ?| ). When hq ^ 0, there exists no pairs of purely imaginary 
eigenvalues A with negative Krein signature, which may bifurcate to complex unstable eigenvalues A 
away from the anti- continuum limit e ^ 0. 

The splitting of the zero eigenvalue of Ti^^\ which defines the stability- instability conclusion of 
Theorem 13.61 may occur in different powers of e as e — > 0. The power of e, where it happens, 
depends on the set S, which classifies the family of the discrete solitons 0„, n g Z. For the sets 
5*1 and S2, which are defined by (|2.13|l and (|2.14|l . we show that the generic splitting of the zero 
eigenvalue occurs in the first and second orders of e, respectively. These results are reported in the 
next two sections. 



4 Bifurcations of the discrete solitons in the set Si 

Here we study the set 5*1 with the explicit perturbation series expansions. These methods illustrate 
the general results of Theorem 13.61 and give asymptotic aproximations for stable and unstable 
eigenvalues of the linearized stability problem H3.2|l . We compare the asymptotic and numerical 
approximations in the simplest cases N — 2 and = 3. 

By Lemma solution of the difference equations H2.8(l is defined by the power series (|2.9(l . where 
(jjn'' is given by 12.5|l with 0„ = {0, tt} for all n G S" and 0^^^ solves the inhomogeneous problem: 

(l-34°)2)4^)=0i"|,+0i^, neZ. (4.1) 

For the set Si, defined by H2.13|) . the system (|4.1|l has the unique solution: 

(/.W = -i(cos((?„_i-0„)+cos(0„+i-0„))e''", 2<n<iV-l, 

= -icos(02-^i)e^'S (Aj^) =-icos(0w-0^_i)e'^~, 

(/.W = e''\ (4-2) 

while all other elements of (jjn'^ are empty. The symmetric matrix Ti. is defined by the power series 
lESJ, where n'-°^ is given by (|3.9I) and Ti^^^ consists of blocks: 
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while all other blocks of Til^Jn are empty. The semi-simple zero eigenvalue of the problem Ti.ip = 7<p 
is split as e > according to the perturbation series expansion: 

(^ = y(") +eV'^'^ +0(e2), 7 = e7i + 0(e2). (4.4) 

Let 7 = be a semi-simple eigenvalue of Ti'-*'-' with N linearly independent eigenvectors f„, n G S*. 
Recalling that sin0„ — and cos0„ = ±1 for all n Cz S, we normalize f„ by the only non-zero block 
(0,cos6'„)"^ at the n-th position, for the sake of convenience. The zero-order term (f^^^^ takes the 
form: 

<^(°) = ^ c„f„, (4.5) 

where c„ G C, ?i G S" are coefhcients of the linear superposition. The first-order term ip^^") is found 
from the inhomogeneous system: 

H(°V^'^ =7i^^"^-W^'V^°^- (4.6) 

Projecting the system (I4.6|l onto the kernel of Ti^^\ we find that the first-order correction 71 is 
defined by the reduced eigenvalue problem: 

Xic = 7ic, (4.7) 
where c = (ci, cn)'^ and Mi is a tri-diagonal A^-by-A^ matrix, given by 

(>^l)m,n= (fm.W^'^fn) , 1 < n, m < iV, (4.8) 

or explicitly, based on the first-order solution H4.2|l and (|4.3|l : 

(Mi)„,„ = cos(6l„+i - On) + cos(6i„_i - On), l<n<N, 
(Mi)„^„+i = (Mi)„+i,„ = - cos(6l„+i - On), l<n<N, 

(Mi)ia = cos(02 - Oi), (Mi)^v,jv - cos{9n - 9n-i). (4.9) 

Similarly, the multiple zero eigenvalue of the problem JHip = AV' is split as e > according to the 
perturbation series expansion: 

^ ^ ,/,(0) + + e,/,(2) + o(e^/i), A = x/iAi + eAa + O(eV^). (4.10) 

Let A = be a multiple eigenvalue of JTiS^^ with N linearly independent eigenvectors f„, n G S" and 
N linearly independent generalized eigenvectors gn, n G S . The eigenector g„ has the only non-zero 
block (cos6'„,0)"^ at the n-th position. The zero-order term is given by H4.5|) as "0'-*''' — while 
the first-order term x/j^^"^ is given by 

'^^'^^yE^-S- (4-11) 

nes 

(2) 

The second-order term tjj^ ' is found from the inhomogeneous system: 

= Ai0(i) + X2iP^"^ - V'"^ • (4.12) 

Projecting the system (|4.12l) onto the kernel of JT-L^^\ we find that the first-order correction Ai is 
defined by the reduced eigenvalue problem: 

2MiC = \\c, (4.13) 
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where A^i is given in H4.8|l . This resuh is in agreement with the leading-order behavior of 
Lemma [3.11 The matrix Aii has the same structure as in the perturbation theory of continuous 
multi-pulse solitons Therefore, the number of positive and negative eigenvalues of A^i is 

defined by the following lemma. 

Lemma 4.1 Let hq, zq, and po be the numbers of negative, zero and positive terms of a„ = 
cos(6'„+i — 9n), 1 < n < N ~ 1, such that uq + zq + pq — N — 1. The matrix M^i, defined by 
14-^^ , has exactly hq negative eigenvalues, zq + 1 zero eigenvalues, and po positive eigenvalues. 

Proof. See Lemma 5.4 and Appendix C of for the proof. ■ 

When Zq — 0, the zero eigenvalue of Aii with the eigenvector (1, 1, 1)-^ is unique. Since all 
On — {0,7r}, n £ S, then all a„ 7^ 0, 1 < j < TV — 1, such that zq — and the splitting of the 
semi-simple zero eigenvalue of 7i^°^ is generic in the first-order of e for the set Si. By Lemma [4. II 
stability and instability of the discrete solitons in the set 6*1 are defined in terms of the number 
Uq of TT-differences in — 0„ for 1 < n < — 1. This result is in agreement with Lemma |3.3I 
and CoroUarv 13.51 for the family Si. Thus, Theorem 13. 61 for the set 5*1 is verified with explicit 
perturbation series results. 

We illustrate the stability results with two elementary examples of the discrete solitons in the set 
5*1: N = 2 and N = 3. In the case N — 2, the discrete two-pulse solitons consist of the Page mode 
(a) and the twisted mode (b) as follows: 

(a) 61=62= 0, (b) 61 =0, 62 = TT. (4.14) 

The eigenvalues of matrix Aii are given explicitly as 71 = and 72 = 2 cos(02 — ^'i)- Therefore, the 
Page mode (a) has one real unstable eigenvalue A ~ 2^/e in the stability problem 1)3. 7|l for small 
e > 0, while the twisted mode (b) has no unstable eigenvalues but a simple pair of purely imaginary 
eigenvalues A ~ ±2iy^ with negative Krein signature. The latter pair may bifurcate to the complex 
plane as a result of the Hamiltonian Hopf bifurcation. 

These results are illustrated in Figures Q and |21 in agreement with numerical computations of 
the full problems H2.8|l and 1)3.2(1 . Fig. ^ shows the Page mode, while Fig. |21 corresponds to the 
twisted mode. The top subplots of each figure show the mode profiles (left) and the spectral plane 
A = Xr + iXi of the linear eigenvalue problem (right) for e = 0.15. The bottom subplots indicate the 
corresponding real (for the Page mode) and imaginary (for the twisted mode) eigenvalues from the 
theory (dashed line) versus the full numerical result (solid line) . We find the agreement between the 
theory and the numerical computation to be excellent in the case of the Page mode (Fig. QJ. For 
the twisted mode (Fig. |2Jl, the agreement is within the 5%-error for e < 0.0258. For larger values 
of e, the difference between the theory and numerics grows. The imaginary eigenvalues collide at 
e « 0.146 with the band edge of the continuous spectrum, such that the real part A^ becomes 
non-zero for e > 0.146. 

In the case iV = 3, the discrete three-pulse solitons consist of the three modes as follows: 

(a) 61=62 = 63 = 0, (b) 61 = 62 = 0, 63 = TT, (c) 61 = 0,62 = TT, 63 = 0. (4.15) 
The eigenvalues of matrix Mi are given explicitly as 71 = and 

72,3 = COS(02 - ^1) + COS(03 ~ 62) ± ^^08^62 -61)- COs(02 - ^l) COs(03 -62)+ 008^(63 -62). 
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Figure 1: The top panel shows the spatial profile of the Page mode (left) and the corresponding 
spectral plane of the linear stability problem (right) for e = 0.15. The bottom subplot shows the 
continuation of the branch from e = to e = 0.15 and real positive eigenvalue theoretically (dashed 
line) and numerically (solid line). 




Figure 2: The top panel shows the twisted mode and the spectral plane for e = 0.15. The bottom 
subplot shows the imaginary and real parts of the eigenvalue with negative Krein signature, which 
bifurcates to the complex plane at e w 0.146. 



12 



Figure 3: Same as Fig. n]for the mode (a) with three excited sites in phase. 
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Figure 4: Same as Fig. El for the mode (b) with three excited sites, where the left and middle sites 
are in phase and the right tt is out of phase. 




Figure 5: Same as Fig. |21for the mode (c) with three excited sites, where adjacent sites are out of 
phase with each other. 
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The mode (a) has two real unstable eigenvalues A w \/6e and A ~ \plk in the stability problem 
ijXTjl for small e > 0. The mode (b) has one real unstable eigenvalue A \/2V3e and a simple 
pair of purely imaginary eigenvalues A « ±i with negative Krein signature. This pair may 

bifurcate to the complex plane as a result of the Hamiltonian Hopf bifurcation. The mode (c) has 
no unstable eigenvalues but two pairs of purely imaginary eigenvalues A ~ ±i\/6e and A ~ ±i\/2e 
with negative Krein signature. The two pairs may bifurcate to the complex plane as a result of the 
two successive Hamiltonian Hopf bifurcations. 

Figures OHSl summarize the results for the three modes (a)-(c), given in H4.15|l . Fig. O corresponds 
to the mode (a), where two real positive eigenvalues give rise to instability for any e 7^ 0. The error 
between theoretical and numerical results is within 5% for e < 0.15 for one real eigenvalue and 
for e < 0.0865 for the other eigenvalue. Similar results are observed on Fig. 0]for the mode (b), 
where the real positive eigenvalue and a pair of imaginary eigenvalues with negative Krein signature 
are generated for e > 0. The imaginary eigenvalue collides with the band edge of the continuous 
spectrum at e « 0.169, which results in the Hamiltonian Hopf bifurcation. Finally, Fig. [Slshows the 
mode (c), where two pairs of imaginary eigenvalues with negative Krein signature exist for e > 0. 
The first Hamiltonian Hopf bifurcation occurs for e w 0.108, while the second one occurs for much 
larger values of e w 0.223, which is beyond the scale of Fig. [S] 



5 Bifurcations of the discrete solitons in the set 5*2 

Here we study the set S2 with the revised perturbation series expansions. The solution is defined 
by the power series H2.9|l . where the zero-order term is given by l|2.5|l with 6n = {0, tt} for all 
n € 5 and the first-order term tpn'' is given by 

^(1) ^ e'''"+i -I- e*^"-i , n^2m, 1 < m < iV - 1, 

4'^=e^'\ 4l^=e'«--^ (5.1) 

while all other elements of (pn'' are empty. The second-order term t/)!^-' solves the inhomogeneous 
problem: 

(1 - 3<^lo)^)0(^) = ^i'l, + 0W , + 3<^(i)^0(°), (5.2) 

with the unique solution: 

= _i (cOs(0„+2 - On) + C0S((?„_2 - ^n) + 2) 6*^" , 71 = 2m - 1 , 2 < m < iV - 1 , 

= -^(cos(^^3-ei) + 2)e^«\ =-i(cos(027V-i-e2W-3) + 2)e'^^"-\ 

0L'l = e'^\ 4'j^+i=e^'— , (5.3) 

(2) 

while all other elements of are empty. The symmetric matrix Ti. is defined by the power series 
H3.8|l . where the zero-order term Ti*^"^ is given by (|3.9|) and the first-order term Ti^^^ is given by 
14.3|l . where = 0, ti G Z. The second-order term Ti^^^ has the structure: 

<) = -20(°V(^) [ n^2m-l,l<m<N (5.4) 
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and 

1 



^'■ni - -'/'i'^' (I ? 1 , n = 2m, < m < iV, (5.5) 



while all other blocks of Ti.n,m are empty. Similarly to the previous section, the semi-simple zero 
eigenvalue of the problem H(p — 7<p is split as e > according to the modified perturbation series 
expansion: 

V = v("'+e¥p(i)+e¥'(2)+0(e3), ^ = e^^^ + Qie"^), (5.6) 
where the zero-order term y^"-* is given by (|4.5|) and the first-order term tp^^^ has the form: 

<p''^ - ^c„(5+f„+5_f„), (5.7) 

neS 

where S± are shift operators of the non-zero 2-block of i„ up and down. The second-order term 
y*^^^ is found from the inhomogeneous system: 

= 72V'"' - T^^^V^ - ?^(2V^°^- (5.8) 
Projecting the system (|5.8|l onto the kernel of Ti^°\ we find the reduced eigenvalue problem: 

M2C = 72C, (5.9) 
where c = (ci, C3, C2n-i)'^ and A^2 is the tri-diagonal N-hy-N matrix, given by 

(■^2)„,„ = (f2™-l,H('^f2„-l) + (f2™-l,7i(^'(5+ +5_)f2„-l) , (5.10) 

for 1 < n, m < A^, or explicitly, based on the first-order and second-order solutions H5.ll) and (|5.3() : 

(M2)n,n = C0S(6I2„+1 - 02n-l) + COs(6'2„-3 " 6'2n-l), 1 < n < TV, 

(M2)«,n+1 = (Af2)„+1,„ = - C0S(6'2„+1 - 6'2„-l), 1 < n < iV, 

(M2)l4 = C0S(6I3 - Oi), iM2)NM = COs(6'2Ar_l - 6'2Ar_3)- (5.11) 

Similarly, the multiple zero eigenvalue of the problem JTixlj — Xip is split as e > according to the 
modified perturbation series expansion: 

= -0'°' + ef/'^'^ + eV^'^ + 0(e3), A = eAi + e^Aj + 0(e3), (5.12) 

where the zero-order term ip^^"^ = y'"-' is given by 14.5|l and the first-order term ■0^^^ has the form: 

The second-order term tp^^"^ is found from the inhomogeneous system: 

jn'^'^^P^^^ = AitA(^) + A2^('') - JH^'^'iP^'^ - ^^(2 V^"^ • (5.14) 
Projecting the system (|5.14|) onto the kernel of Jl-L^^\ we find the reduced eigenvalue problem: 

27W2C = A?c, (5.15) 

15 
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Figure 8: Same as Fig. Obut for the mode (a) of the set 82- 
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in accordance with Lemma 13.11 Since the matrix A4 2 has exactly the structure of the matrix A4 1 , 
described in Lemma 14.11 we conclude that the stability and instability of the discrete solitons in 
the set 5*2 is defined in terms of the number no of 7r-differences in 02n+i — ^2n-i, I < n < N ~ 1, 
in accordance with Lemma and Corollarv l3.5l Thus, Thcorcm l3.6l for the set S2 is verified with 
explicit perturbation series results. 

We summarize that bifurcations and stability of the discrete solitons in the set S2 is exactly equiv- 
alent to those in the set Si, but the splitting of all zero eigenvalues occurs in the order of e^, rather 
than in the order of e. These results for the set 6*2 with N = 2 and = 3 are shown on Figures 
I ^IUI in full analogy with those for the set Si. The corresponding asymptotic approximations of 
eigenvalues can be "translated" from those of the previous section by substituting -/i e. Fig. 
shows the Page mode where the agreement with the theory is excellent for e < 0.2. Fig. shows 
the twisted mode with very good agreement for e < 0.415 and the Hamiltonian Hopf bifurcation at 
e « 0.431. The only difference from the twisted mode of Fig. |21is that the imaginary eigenvalue of 
negative Krein signature collides with the imaginary eigenvalue of positive Krein signature, rather 
than with the band edge of the continuous spectrum. Figs. and 1101 show the modes (a), (b), 

and (c), respectively, of the three excited sites. Again, the Hamiltonian Hopf bifurcations occur 
when the imaginary eigenvalues of negative Krein signature collide with the imaginary eigenvalues 
of positive Krein signature. For the mode (b), the bifurcation occurs at e « 0.328 (see Fig. EJ. For 
the mode (c), two bifurcations occur at e ^ 0.375 and e « 0.548 (see Fig. [TU)l . 

6 Summary 

We have studied stability of discrete solitons in the one-dimensional NLS lattice (|1.1|) with d = 1. 
We have rigorously proved the numerical conjecture that the discrete solitons with anti-phase excited 
nodes are stable near the anti-continuum limit, while all other discrete solitons are linearly unstable 
with real positive eigenvalues in the stability problem. Additionally, we gave a precise count of the 
real eigenvalues and pairs of imaginary eigenvalues with negative Krein signature. These results 
are not affected if the excited nodes are separated by an arbitrary sequence of empty nodes. We 
studied two particular sets of discrete solitons with explicit perturbation series expansions and 
numerical approximations and found very good agreement between the asymptotic and numerical 
computations. 

Stability and instability results remain invariant if the discrete solitons are excited in the two- 
dimensional NLS lattice Hl.l|l with d = 2, such that the set S is an open discrete contour on the 
plane. Similar perturbation series expansions for the sets Si and ^2 in the two-dimensional NLS 
lattice can be developed and the same matrices Aii and AI2 define stability and instability of these 
discrete solitons. 

In the second forthcoming paper of this series, we shall consider a closed discrete contour on the 
plane for the set S. Such sets may support both discrete solitons and discrete vortices with a non- 
zero topological charge. Continuation of the limiting solutions from e = 0toe^0isa non-trivial 
problem if the amplitudes (/>„ are complex- valued. We shall study persistence, multiplicity and 
stability of such continuations with the methods of Lyapunov-Schmidt reductions. 
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